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ABSTRACT 

A simple model of the lower corona which allows for a possible 
difference in the electron and proton temperatures is analyzed. With 
the introduction of a phenomenological heating term, temperature and 
density profiles are calculated for several different cases. It is 
found that, under certain circumstances, the electron and proton 
temperatures may differ significantly. 
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I . INTRODUCTION 


It is generally believed that the corona is heated by a flux of 
non-thermal energy which originates in the photosphere and propagates 
through the chromosphere, transition region, and possibly part of the 
corona. D'Angelo (1969) suggests that if this non-thermal flux is in 
the form of ion acoustic waves, Landau damping of these waves might 
contribute significantly to the heating of the corona. D'Angelo also 
noted that if this were the case, preferential heating of the protons 
could cause the proton temperature to be higher than the electron 
temperature . 

There are several interdependent aspects of the problem of coronal 
heating. These include the propagation of non-thermal energy, the 
dissipation of this energy, and the effect of this dissipation on the 
temperature and density structure of the corona. In this article, we 
shall be concerned only with the last topic, and we shall consider 
explicitly the possibility that a substantial fraction of the dissipated 
energy preferentially heats the ions. 

II. ANALYSIS 

We adopt a simple planar model of the corona in which all quantities 
are uniform in the horizontal coordinates and in which the gravitational 
acceleration is independent of the vertical coordinate (z). The composition 
is taken to be that of fully ionized hydrogen gas. Radiation losses and 
magnetic fields are neglected. Since we are not treating the problem of 
dissipation at this time, a phenomenological heating term is introduced. 

We suppose that all of the dissipation occurs in a narrow region, of 
width of order 10,000 kilometers, centered at a height of 15,000 kilometers 
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above the base of our model, which we take to be the level at which 
T = 10^ °K . The particular form adopted for the heating function 

^ , ” 3 -lx 

Q (erg cm s ) is 



z, £ z <. z„ + z 
1 1 w 


( 1 ) 


where = 10,000 km. Q = 0 outside the range given in equation (1) . 

We further assume that a fraction f of the total energy input goes to 
heating the protons, so that 


Q_ = (1 -f) Q, 


Q - fQ 

+ 


( 2 ) 


where Q (z), Q (z) are the electron and proton heating functions, 
respectively. 

The relevant equations will be taken from Braginskii (1965). The 


equations of hydrostatic balance are 


— + en E - mng = 0 

dz + + + 


-dp_ 

dz 


-enE-mng^o, 


(3) 

(4) 


where the symbols have their usual meaning. The subscripts + and 

refer to protons and electrons, respectively. We may set n + n_ = n 

in the above equations (but not in Poisson’s equation). Then, neglecting 

m in comparison with m and summing the above equations, we obtain 

+ 


h (»k[T + + T j) + nm + 


g = 0 . 


(5) 


Two heat equations may be expressed as 


d 

dz 


(h _t)+Q - | v nk (T - T 

\ + dz / + 2E + 


T ) = 0 


( 6 ) 


and 
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d 

dz 


(■- +) 


+ Q 


+ ^ v £ nk (T + - T_) ^ 0 . 


(7) 


The thermal conduction coefficients h , h_ and the heat-exchange 


coefficient v are given by 

£j 


-4.73 -1 5/2 -6.12 -1 _ 

k = 10 X T , k = 10 X T 

+ + 


5/2 


-2.41 , -1 -3/2 

V = 10 X n T_ , . (8) 

E 


where X (the Coulomb logarithm) can be expressed with sufficient accuracy 
as 

\ = 15.9 - 1.15 log n + 2.3 log T . (9) 


For the ranges of density and temperature of interest, we may with 
sufficient approximation adopt X = 20. Then the coefficients in 
equation (8) may be expressed as 


7 5/2 7 5/2 2 -1 3/2 

= ? K _T_ - %=2 K + T + ’ V E = 3 Ak n T - 


( 10 ) 


-6.57 7.96 -16.79 

where K = 10 , K = 10 and A = 10 

+ 

With these substitutions, equations (6) and (7) become 

d T - = - K _1 Q - AK -1 n 2 T" 3/2 (T - T ) , 
dz + 


(ID 


dT + =-k 1 Q + ak” 1 n 2 T 3/2 (T - T ) , 

■ + + + - + 

dz 


( 12 ) 


from which we see that 

4 ( k - t - 7/2 + Vl /2 ) - -Q • 

dz 

Equation (5) now becomes 

I — = - (T + T )“ 1 (C + - [T + T ]) 

n dz + dz - + 

where 


(13) 


( 14 ) 
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( 15 ) 


„ , -1 - _-3 . 48 

C = m gk pa 10 
+ 

III. RESULTS 

Our aim is to calculate n(z), T_(z) and T (z) for specified 

forms of Q (z) and Q (z), and for given boundary conditions. One 

+ 

of the five required boundar}^ conditions is as follows, 

n =5 n at z = 0 , (16) 


Since equation (13) is immediately integrable, it is convenient to 
express two of the conditions as follows, 


K T 


7/2 


K T 
+ + 


7/2 


- (K_+ 


V 


7/2 


at z = 0 


(17) 


and 


(K 

dz 


,7/2 


7/2 

+ K T* ) 
+ + 



( 18 ) 


Equation (18) implies that the entire energy input flows back down as 
a heat flux. The remaining conditions are taken to be the following: 


T - T_ -♦ 0 , z - <» , 

“f r 


T - T -* 0 , z , 

+ 


(19) 

( 20 ) 


In selecting the boundary conditions at the base of the model, we 
have considered estimates for the product nT in the transition region. 

14 

Athay (1965) estimated this quantity to be in the range of (3. 7-6.0) x 10 ; 

14 

later, Athay (1969) made the estimate 6 x 10 ; Moore and Fung (1972) 

14 

adopt the range (5-15,1) x 10 . From these estimates, we have chosen 

6 

the values log n = 8.50, 8.78, 9.00 at T = 10 . We also consider 
o 

-2 -i 5 

two values of the downward heat flux F (erg cm s ): 6 x 10 
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(Dupree and Goldberg 1967), and 10 (Moore 1972). 

The results of our calculations are displayed in Figures 1 and 2. 

We see from Figure 1 that an increase in leads to a similar increase 

in n(z) , so that the profile is insensitive to n . An increase of 

o 

n hardly affects T (z), but it leads to a reduction of T (z) since 
o + 

it tightens the coupling between the two species. (Note, however, that 

the asymptotic value of T , as z -» ® , is not significantly affected 

+ 

by n since T -+ T in this limit.) 
o + - 

Figure 1 also shows that, if n Q and F are fixed, an increase 

in f leads to an increase in T + in the heating region and to a 

slight reduction in n in this region. The decrease in density tends 

to offset the increase in pressure due to the increase in proton temperature. 

The electron temperature is not significantly effected by a change in 

f. We see from Figure 2 that, for fixed F and f, an increase in n^ 

leads to an increase of n(z) and a reduction in T (z) , as noted 

+ 

above, but this diagram shows that there is a slight increase in T_(z) 

in the heating region. A comparison of Figures 1 and 2 shows that an 

increase in F, for fixed n , f, leads to increases in n(z), T (z) 

o 

and T (z). 

+ 

8.5 

It is notable that, if n is as small as 10 , T may exceed 

o + 

T , in the heating region, by up to 70%. 

The above calculations show that, if the non-thermal energy goes 
primarily into heating the proton gas, the proton temperature may be 
substantially higher than the electron temperature in the heating region 
of the corona. Hence one would need to be cautious in inferring the 
coronal structure from observational data, since some observations 
(such as radio-frequency brightness and coronal line strengths) depend 
upon electron temperature, other observations (such as coronal -line widths) 
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depend upon the ion temperature, and some observations (such as the 
variation of density with radius) depend upon a combination of the 
temperatures . 

In order to pursue this question, it is essential to investigate 
iti more detail the heating of the solar corona, to determine in particular 
the allocation of energy between the electron gas and the proton gas. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2 . 


Temperature and density profiles for a total energy input of 
6 -2 —I 

1.0 x 10 erg cm sec . Case I corresponds to an initial 
9 -3 

density of 10 cm and f = 1. Case II also corresponds to 

9 —3 

an initial density of 10 cm but with f = 0.5. Case III 

8 . 78 -3 

corresponds to an initial density of 10 cm and f = 1. 

The T profiles are indistinguishable on this scale. 

Temperature and density profiles for a total energy input of 
5 -2 -1 

6 x 10 erg cm sec and f = 1. Case I corresponds to an 

9 -3 

initial density of 10 cm . Case II corresponds to an 

8 78 “3 

initial density of 10 * cm (The T_ curves for these 

first two cases are indistinguishable.) Case III corresponds 

8 5 —3 

to an initial density of 10 cm 
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